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ABSTRACT 

The tidal disruption of a star by a supermassive black hole provides us with a rare glimpse 
of these otherwise dormant beasts. It has long been predicted that the disruption will be ac- 
companied by a thermal 'flare', powered by the accretion of bound stellar debris. Several 
candidate disruptions have been discovered in this manner at optical, UV and X-ray wave- 
lengths. Here we explore the observational consequences if a modest fraction of the accretion 
power is channeled into an ultra-relativistic outflow. We show that a relativistic jet decelerates 
due to its interaction with the interstellar medium at sub-parsec distances from the black hole. 
Synchrotron radiation from electrons accelerated by the reverse shock powers a bright radio- 
infrared transient that peaks on a timescale ~ 1 yr after disruption. Emission from the forward 
shock may be detectable for several years after the peak. Deep radio follow-up observations of 
tidal disruption candidates at late times can test for the presence of relativistic ejecta. Upcom- 
ing radio transient surveys may independently discover tens to hundreds of tidal disruptions 
per year, complimenting searches at other wavelengths. Non-thermal emission from tidal dis- 
ruption probes the physics of jet formation under relatively clean conditions, in which the flow 
parameters are independently constrained. 

Key words: black hole physics - galaxies; nuclei. 
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1 INTRODUCTION 

Supermassive black holes (SMBHs) are most easily studied when 
they accrete at high rates for extended periods of time and power 
active galactic nuclei (AGN). However, the majority of galactic nu- 
clei are relatively quiet. With the exception of Sgr A* and a hand- 
ful of nearby low luminosity AGN, underfed SMBHs are difficult 
to study. Even in bright AGN, obtaining a detailed understanding 
of the accretion process is hindered by our incomplete knowledge 
of the 'boundary conditions' of the flow, such as the accretion rate 
and the magnetic field strength. 

A rare glimpse into the properties of normally quiescent 
SMBHs is afforded when a star passes sufficiently close that it is 
torn apart by the SMBH's tidal gravitational field (Hills 1975; Rees 
1988; Goodman & Lee 1989). Analytic estimates and numerical 
calculations show that the process of disruption leaves a significant 
fraction of the shredded star gravitationally bound to the black hole 
(e.g. Rees 1988; Ayal et al. 2000). The accretion of this stellar de- 
bris has long been predicted to power a thermal 'flare' at optical, 
UV, and X-ray wavelengths that lasts for months to years (e.g. Ul- 
mer 1999; Stubbe & Quataert 2009, 2010). 

Models of the stellar dynamics in galactic nuclei estimate that 
the tidal disruption rate is ~ 10"^^ - lO^^yr"' per galaxy (Magor- 
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rian & Tremaine 1999; Wang & Merritt 1004). Although tidal dis- 
ruptions are rare, and potentially difficult to detect in the bright 
and dust-extincted central regions of galaxies, surveys at X-ray 
(e.g. Donley et al. 2002; Esquej et al. 2008), optical (van Velzen 
et al. 2010) and far UV wavelengths (Gezari et al. 2006, 2008, 
2009) have now detected ~ 10 candidates (see Gezari 2009 for 
a recent review). Future wide-field surveys, particularly at optical 
wavelengths (e.g. the Large Synoptic Survey Telescope), have the 
potential to detect hundreds of tidal disruptions per year (Strubbe 
& Quataert 2009). A detailed census of these events would pro- 
vide a unique probe of SMBH demographics and the shape of the 
gravitational potential in galactic nuclei (e.g. Merrit & Poon 2004). 



In this Letter we consider a different observational signature 
of tidal disruption that occurs if the accretion of stellar debris pow- 
ers a transient relativistic jet (^. We show that if a modest frac- 
tion of the accretion power is used to accelerate material to ultra- 
relativistic speeds, the energy that is released when the jet interacts 
with the interstellar medium (^ may power bright non- thermal 
synchrotron emission at radio-infrared wavelengths that peaks on a 
timescale ~ 1 year after disruption (®. We explore the prospects 
for detecting such radio flares in ^and discuss our results in ^ 
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2 STELLAR DISRUPTION AND ACCRETION 

Disruption occurs when the tidal force apphed by the SMBH over- 
comes the self gravity of the star. This occurs when the radius of 
orbital pericenter i?p is less than a critical distance (Rees 1988) 



(1) 



where Mj, and R^, are the stellar mass and radius, respectively. 
For black holes masses Mbh <: 10^ Mq and solar-type stars with 
radii if* ~ 10" cm, disruption occurs inside the event horizon, 
i.e. Rf ^ Rs ~ '^GM^n/c-, such that no electromagnetic signa- 
ture is expected. For less massive black holes, disruption occurs 
outside Rg. However, the prevalence of very low mass SMBHs 
(^BH ^ 10^ Mq) in galactic nuclei is unclear (e.g. Ferrarese et 
al. 2006; Greene & Ho 2007) and, even if they are abundant, their 
electromagnetic emission will be dim if it is limited to a fraction of 
the Eddington luminosity L^m ^bh- The 'sweet spot' for bright 
tidal disruption, therefore, occurs for Mbh ~ 10'"** Mq. 

Disruption unbinds ~ 1/2 of the stellar debris from the sys- 
tem (e.g. Rees 1988). The rest of the mass is placed onto highly 
eccentric orbits, which return it to the vicinity of the black hole on 
a wide range of timescales (Lacy et al. 1982; Evans & Kochanek 
1989; Laguna et al. 1993). The most tightly bound material returns 
on a timescale (e.g. Strubbe & Quataert 2009) 
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where the last expression is evaluated for a solar-type star (Rj, = 
/?o) with a pericenter distance = 6GMbh/c^ $ Rt- For clarity we 
focus on this case throughout the remainder of this Letter. 

As material returns to the SMBH, it shocks on itself, circular- 
izes, and accretes (e.g. Kochanek 1994). The accretion rate peaks 
at / ~ ^fallback and declines oc r^'^ thereafte£] (Rees 1988), viz. 
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For Mbh ^ 4 x 10' Mq the accretion rate is initially super- 
Eddington, i.e. M > Mgdci = WLem/c^ a 0.2(Mbh/ 10' M0)Mo 
yr"' . Accretion decreases below the Eddington rate on a timescale 



fEdd-0.24(MBH/10'Mof' yr 



(4) 



A ~ 10' Mq black hole with an accretion rate M ~ Mem produces 
thermal emission that peaks in the Far-UV {kT ~ 30 eV), but which 
may also be detectable at optical, near-UV, and X-ray wavelengths. 
Thermal 'flares' with a characteristic duration ~ t^M ~ months are 
thus considered a hallmark signature of tidal disruption; their blue 
optical-UV colors are the primary means by which they are cur- 
rently identified (Gezari et al. 2009). 

2.1 Relativistic Jet Production 

Black hole accretion is often accompanied by powerful, collimated 
relativistic outflows. A comparison of the disk luminosity and jet 
power of blazars, for instance, shows that a substantial fraction 
of the accretion power goes into relativistic jets (e.g. Ghisellini et 
al. 2009). A similar correlation between jet and accretion power 
is observed in elliptical galaxies (Churazov et al. 2002; Allen et 
al. 2006). 

Though common, powerful jets are not present in all accreting 



' Note, however, that corrections to this power-law decline may occur at 
both early (Lodato et al. 2009) and late (Cannizzo et al. 1990) times. 



systems at all times. In stellar mass X-ray binaries, for instance, 
jets are ubiquitous at low accretion rates M < 0.03MEdd (the low- 
hard state), but they are intermittent (or absent entirely) at higher 
accretion rates (the 'thermal/high' and 'very high' states; Fender et 
al. 2004). The typical radio 'loudness' of AGN also increases for 
smaller Eddington ratios A = MIMeaa^ but nevertheless some frac- 
tion of AGN show evidence for jets at all values of A. Furthermore, 
the 'radio loudness dichotomy' (e.g. Ho 2002; Sikora et al. 2007) 
suggests that a 'second parameter' other than M, such as the black 
hole spin (Wilson & Colbert 1995) or the magnetic flux threading 
the disk (Spruit & Uzdensky 2005), may also control the jet power 

The total energy released by the accretion of a solar mass star 
is iiacc ~ 10^3 erg for an assumed accretion efficiency ~ 0. 1. The 
fraction of this energy ej placed into relativistic ejecta is highly un- 
certain and could vary between disruption events, depending e.g. on 
the magnitude of the SMBH spin and its direction with respect to 
the angular momentum of the disk. We leave 6j as a free parameter 
in our calculations, but adopt a fiducial value ej = 0.01 as motivated 
below. The total energy of the jet is thus £j = Cjiiacc ~ 10^'(ej/0.01) 
erg, similar to that of a supernova or gamma-ray burst (GRB). We 
furthermore assume that the jet has an opening angle » 0. 1 and 
a bulk Lorentz factor Fj = 10, values typical of AGN jets. 

We consider two scenarios for producing a relativistic jet fol- 
lowing tidal disruption. The first assumes that a relativistic out- 
flow is present during the early super-Eddington accretion phase 
(M > MEdd; t < fEdd ~ 10^ s; eq. Q), and that during this 
time the jet kinetic luminosity is constant Lj = ^j/fj, where 
tj ~ 'Edd is the jet duration. The jet efficiency can be written 
as fj ~ O.l(LJ/LEdd)(MBH/lO'M0)'/^ such that 6| = 0.01 for 
Lj = 0. ILEddi Mbh ~ 10' Mq. On one hand, this scenario is conser- 
vative because numerical simulations of super-Eddington accretion 
show that the accretion rate reaching the black hole can far exceed 
Mfidd (e.g. Ohsuga et al. 2005), in which case Lj could in principle 
be » Lfidd- On the other hand, whether a relativistic jet accompa- 
nies super-Eddington accretion is observationally uncertain. 

In a second, more conservative scenario we assume that a rel- 
ativistic jet forms only once the accretion rate M oc r^'' (eq. ||3]) 
decreases to < 0.03MEdd = M, as occurs at times t > lOfsdd- When 
the accretion rate becomes this low, the accretion disk is predicted 
to become radiatively inefficient, geometrically thick, and particu- 
larly susceptible to outflows (e.g. Narayan & Yi 1995; Blandford 
& Begelman 1999), a fact supported by observations of X-ray bi- 
naries in the low/hard state and the radio loudness of AGN at low 
Eddington ratios. The total energy of the jet power in this scenario 
is £j ~ 0.lM,c-fj where the jet 'duration' is fj ~ lOJEdd- The jet 
efficiency can be written e, a 0.03(Lj/0.1Mc-)(Mbh/10^Mo)'/\ 



3 JET-ISM INTERACTION 

Jets from AGN are long-lived and propagate to ~kpc-Mpc dis- 
tances before dissipating their bulk energy (powering e.g. giant 
radio lobes). By contrast, any relativistic outflows from tidal dis- 
ruption are transient and hence decelerate via interaction with the 
interstellar medium (ISM) at a much smaller distance < pc. As the 
jet drives a forward shock (FS) into the ISM, a reverse shock (RS) 
propagates through the ejecta, slowing it down ( 33.11 1. As discussed 
below, once the reverse shock crosses the jet, the entire configura- 
tion (swept-up ISM and shocked ejecta) has been slowed to mildly 
relativistic speeds. After this point, the flow relaxes into a quasi- 
spherical, non-relativistic Sedov phase 03.2I see Fig.[T]l. 
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Figure 1. Sketch of the ultra-relativistic (left) and Newtonian (right) stages 
of jet deceleration. Initially the jet propagates radially, forming a for- 
ward/reverse shock (FS/RS) structure as it interacts with the ISM ( j|3.U . 
Once the configuration has slowed to non-relativistic speeds, each of the two 
jet lobes expands quasi-spherically ( ^3.2) . For a typical observer, emission 
peaks at the transition between the relativistic and non-relativistic stages 
due to synchrotron emitting electrons accelerated by the reverse shock 
(shown as grey shading on the right). 



3.1 Relativistic Phase 

The interaction between a relativistic jet and the ISM has been stud- 
ied analytically in two limits. In the first approach, one consid- 
ers a quasi-steady jet that slowly propagates through the ambient 
medium, creating a characteristic jet-cocoon structure (e.g. Begel- 
man & Cioffi 1989). This limit is appropriate for long-lived AGN 
jets. A second approach considers an impulsive ejection of rela- 
tivistic material (a thin shell or 'pancake') that produces a forward- 
reverse shock structure upon interacting with the ISM (Sari & Piran 
1995). This limit applies to GRB jets. 

Neither the steady nor the impulsive limits strictly apply to 
the jets from tidal disruption. We nevertheless adopt the impulsive 
limit because the time required for the jet to appreciably deceler- 
ate ~ ^cr/c (defined below) is > to the jet duration fj. We fur- 
thermore approximate the intrinsically multidimensional problem 
as the one dimensional (ID radial) deceleration of a conical jet. 
This approach is supported by relativistic hydrodynamical simula- 
tions (e.g. Zhang & MacFadyen 2009), which demonstrate that 2D 
efi^ects, such as sideways expansion, become important only in the 
late, non-relativistic phases of the interaction ( 93. 2b . 

As the jet drives a FS into the ISM, a RS propagates back into 
the unshocked ejecta. By assuming pressure balance at the contact 
discontinuity and solving the shock jump conditions, Sari & Piran 
(1995) calculate the bulk Lorentz factor Fsh of the shocked fluid as 
a function of the radius R of the flow. For ISM of constant density 
nisM they show that while the reverse shock crosses the ejecta, the 
shocked fluid decelerates as Fsh(/^) - Fci(i?/i?cr) where 

«„^.r;-,„.,.,„"(^)'"(^)-"'(^)"'e™,s, 

is the radius at which the reverse shock fully crosses the ejecta, and 



F„ = 



is the Lorentz factor at R 
jet energy and fi, = '^^ 



^sh- Here Ei^a = Ejf^ is the 'isotropic' 
(?/2 a; 5 X 10"^^ is the jet beaming fraction. 



We scale the ISM density to ~ 10 cm ^, similar to that inferred 
from X-ray observations of the inner parsec of the Galactic centre 
(e.g. Baganov et al. 2003; Quataert 2004). 

Equations ^ and l|6) show that the jet decelerates to a mildly 
relativistic velocity at R^ ~ 0.1 - 1 pc on a timescale /ci ~ 
^ci7/8crC ~ 1 and 3 years for /j ~ and ~ lOfedd. respectively. 
Note that this result depends very weakly on the precise values of 
6j and «isM- Importantly, most of the bulk kinetic energy of the jet is 
dissipated by the reverse shock, which decelerates the ejecta from 
its initial Lorentz factor Fj ~ 10 to its final value F^ ~ 1 with most 
of the dissipation occuring at radii ~ Dissipation by the for- 
ward shock is sub-dominate because it is only mildly relativisti(|| 
with Fps ~ Fc, ~ i diXR ~ R^y. As we discuss in ^ emission from 
the jet is also dominated by the reverse shock. 

3.2 Non-Relativistic Phase 

After the RS has crossed the jet (R > R,.y), the fluid expands mildly 
relativistically. After a transient phase the outflow along each jet 
head relaxes into a non-relativistic, quasi-spherical expansion. The 
evolution of each 'lobe' at late times can thus be approximated by 
a Sedov-Taylor solution, for which the radial velocity evolves as 
j8nr ~ r^^'^, where a lower case r denotes the radial distance mea- 
sured from the center of each lobe (at radius R ~ R^ from the 
black hole; see Fig.[T]for an illustration). During the non-relativistic 
phase the gas pressure decreases as P oc /3^^ oc r"^, such that adia- 
batic expansion cools the electrons accelerated at the reverse shock 
(denoted by the shaded region in Fig.[T}. As discussed in 34.21 this 
controls the rate at which the RS emission declines after the peak. 

4 EMISSION 

In calculating the synchrotron emission from the forward and re- 
verse shocks, we make several standard assumptions regarding the 
shock microphysics (e.g. Sari et al. 1998), as motivated by the phe- 
nomenology of GRB afterglows. First, we assume that a fraction 
6e ~ 0.1 of the energy dissipated at the shock is used to accel- 
erate relativistic electrons. We assume the electron distribution is 
a power-law with index p ~ 2.5, such that the minimum elec- 
tron Lorentz factor is y,,, = e(,(F,ei - l)(p - 2)m,,/(/j - l)mg x 
600e(,(Frci - 1), where F^i is the relative Lorentz factor between 
the shocked and unshocked fluid. We furthermore assume that the 
magnetic energy density of the shocked fluid is a fraction eg ~ 10"^ 
of the total energy density!^ 

4.1 Peak Emission 

In general the axis of the jet will be directed an angle dobs ~ 1 rad 
with respect to the line of site. Emission is thus weak during the 
initial relativistic stage 03.U because (1) the luminosity is Doppler 
deboosted by a factor oc 6*, where 6 = T^^{1 - j8sh cos (^obs) ' ~ 
2/Fsh <K 1 and (2) only a modest fraction of the jet has interacted 
at times t <s: t^ (R Ra)- Emission peaks when the reverse shock 
passes entirely through the ejecta {t ~ fpea^ ~ ?„; R ~ Rcr) because 



^ In GRB jets the opposite holds true; the forward shock is relativistic and 
the reverse shock is mildly relativistic. This is because the duration of the 
GRB fGRB « (j. 

^ If the jet is strongly magnetized at i? ~ R^, then both the particle accel- 
eration (e.g. Sironi & Spitkovsky 2009) and dissipation at the reverse shock 
(Mimica et al. 2009) may be strongly suppressed. However, since decel- 
eration occurs at large radii Rcr ~ 10^"^i?s, it is reasonable to assume that 
jet has fully accelerated, in which case the residual magnetization may be 
sufficiently low for a strong shock and non-thermal particle acceleration. 
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a large fraction of the jet energy is dissipated around this time and 
the expansion has become mildly relativistic {6 ~ 1). 

We now estimate the synchrotron emission at peak brightness. 
Motivated by the above discussion, in what follows we assume that 
ySsh = 0.6 (Fsh = 1.25) at f ~ t^.j^. On timescales fp^ak ~ ^ci^ the 
relative Lorentz factor between the shocked and unshocked fluid is 
TreLRs ~ rjr,,(l - rj/2 and r,ei,Fs = for the RS and 

FS, respectively. The minimum Lorentz factor of the accelerated 
electrons is 7n,,FS ^ 20(6,/0.1) and 7„,.rs ^ 300(Q/0.1)(rj/10). 
Synchrotron flux peaks at the frequency 



eBjj . 0.l(^)-(^)' (^)' GHz (FS) 
2nm,c I 25(^f(jg,)"'(^)"^(Iif GHz (RS) 



where B = yfSne^ and eth - ('^^a + ^Wsh - 1)"',,«ismc^ - 
ImpiiisMC^ are the magnetic field strength and energy density of the 
shocked fluid, respectively (e.g. Sari & Piran 1995). 

The synchrotron luminosity at Vn, is L,ot - (4/3)A^eCrrC7^L's, 
where A^^ is the number of radiating electrons and Ub = B^/Sn. 
At ?peak ~ ?cr the numbcr of ISM electrons swept up by the FS 
A^e Ps = 2n6jRl^nisyi/3 ~ 10^^ is comparable to the number in the 
shocked jet iVe.RS = ^j/WpC^Fj ~ 7 x W^{eJQm){TJW)-K The 
total emission from the RS thus exceeds that from the FS by a factor 
~ (rm.Rs/rm.Fs)^ > 10^ (Fj / 10)^ Although naivcly the FS might Still 
appear to dominate at low frequencies, self-absorption suppresses 
the emission at v < 1 GHz. The contribution of the FS is modest at 
the peak emission. The synchrotron luminosity Ltot of the RS is 



A„, « 6 X 10*''(ej/0.01)(e,/0.1)2(6s/0.01)(Fj/10) ergs" 



(8) 



The peak flux at v « Vn, for an event at distance D is thus given by 
fv„ - Ltot/'i^D\^ or 



F ,/_^U_fs_V'V_;^isM_y'- 

I1O-2/I1O-2/ llOcm-V 



10/ \Gpc/ 



mJy, (9) 



Note that if the jet is modestly energetic (ej ~ 10"'' - 10"^), the 
predicted flux ~ 0. 1 - 1 mJy is well within the sensitivity of cuiTent 
radio telescopes at ~Gpc distances (redshift z ~ 0.2) similar to 
those of candidate disruption events. 

Figure |2] shows the spectrum of the RS emission at f a ?p(;ak ~ 
t„. In addition to the synchrotron peak frequency at v^, « 30 GHz 
(eq. I?]), the two other breaks denote the self-absorption and cool- 
ing frequencies (at lower and higher frequencies, respectively). Self 
absorption occurs below the frequency at which the optically-thin 
synchrotron luminosity equals the Rayleigh- Jeans luminosity of a 
plasma with temperature = yc^RsnigC^ /3. The cooling frequency 
is the characteristic synchrotron frequency of an electron with a 
cooling timescale equal to the expansion timescale ~ f^r. 

In Figure |2] we plot for comparison an estimate of the con- 
temporaneous emission from the accretion disk (assuming an ex- 
tended thin disk). At / ~ fpe;,^ ~ 1 year the disk luminosity is 
Li ~ 0.lMc\,^ ~ 10*^ erg s"' for Mbh ~ 10' Mq, i.e. a fac- 
tor ~ 10^ brighter than the that of the RS for fiducial parameters 
(eq. (8)). However, because the disk emission peaks in the Far-UV, 
emission from the jet dominates at radio-IR wavelengths. Although 
emission from the disk appears to dominate that from the jet at 
optical-UV wavelengths, the opposite may be possible for more 
realistic disk models. We note that the light curves of a few tidal 
disruption candidates 'flatten' on a timescale f ~ 1 year (e.g. Dl-9, 
D3-13; Gezari et al. 2008, 2009). If this flattening is due to the con- 
tamination of emission from the jet, we predict that the optical-UV 
spectrum should simultaneously flatten to vFy ~ const (see also 
Wong et al. 2007). 
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Figure 2. Predicted synchrotron spectram (solid line) at the time of peak 
emission (peak ~ 'cr ~ 1 year from the tidal disruption of a solar mass star by 
a 10' Mq black hole at a distance D = I Gpc. From low to high frequencies, 
the three breaks correspond to the self absorption, characteristic and cooling 
frequency, respectively. The dashed line shows an estimate of the themial 
emission from the disk at a similar epoch. 



4.2 Light Curve Rise and Decay 

In this section we estimate how the light curve rises prior to, and 
declines following, the epoch of peak emission {t ~ tp^ak ~ fcr)- 
Prior to the peak t <c t^^, the RS is still crossing the jet, and the 
Lorentz factor of the shocked fluid decreases as Fsh 0^ Rr^^~. The 
characteristic synchrotron frequency and the radiated power scale 
as Vn, oc 6By^ and L,o, oc S^N^y^B-, respectively. For a typical 
observer angle 9a^,^ ~ 1 rad, the Doppler factor evolves as (5 oc 
T^2^' oc while the time measured by the observer is tohs. ~ 
(R/pshC)(i - Ah COS Sobs) ~ R/2c oc R. The energy density of the 
shocked fluid oc B^ scales as oc F^i^ oz R \ the minimum electron 
Lorentz factor y^ oc Fj/Fsh oc if''^, and the number of electrons 
swept up by the reverse shock evolves as Ne oc J(ySj - Psh)dR ~ 
^.sh/Fji, oc R^^. Combining these results, we find that Vy„ oc R oc f^i,,, 
R* oc f*^^, and F„, oc /^tot/^m ^ ^obs' '^his Steep predicted rise 
in the radio flux suggests that observations at times <c fp^ak ~ 1 yr 
might be too early to detect the radio emission. 

Once the RS crosses the jet, only the FS continues to accel- 
erate electrons. Nevertheless, for some period after the peak, the 
emission is still dominated by the electrons originally accelerated 
by the reverse shock when R ~ R^^, which adiabatically cool as 
the outflow expands. Once the fluid transitions to a non-relativistic 
Sedov-Taylor expansion, its radius increases as r oc and its pres- 
sure decreases as P oc r^^; the temperature of the relativistic elec- 
tron 'gas' Tg oc y^ thus decreases oc oc r"^!^ . The magnetic 
field generated by the shock decreases B' oc P oc r^^. Combining 



these results, one finds that v„. 



obs ' and 'oh. 



3/5 



This relatively steep decay, combined with the steep rise predicted 
above, suggests that the duration of the peak St is « t^cak- Hy- 
drodynamic simulations of the trans-relativistic stage are, however, 
required to more precisely quantify the light curve shape. 

Although RS electrons dominate the emission immediately af- 
ter fpeak, synchrotron radiation from the 'fresh' electrons acceler- 
ated by the FS will take over at later times. A useful analogy in this 
case are the radio afterglows of nearby GRBs. For GRB 030329, the 
redshift z - 0.17 and total energy in relativistic ejecta ~ 3 x 10^" 
ergs are similar to the characteristic values for tidal disruption. The 
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2.3 GHz flux of 030329 was Fy ~ mjy at f ~ 1 year after the burst, 
with the light curve decaying relatively gradually (Fy oc r'; van 
der Horst et al. 2008). 

5 DETECTION PROSPECTS 

One method to detect transient radio emission from tidal disruption 
is via follow-up observations of candidate events detected at other 
wavelengths (e.g. optical, UV and X-rays). The only relevant radio 
follow-up to date (of which we are aware) is of the SDSS Stripe 82 
candidate TDE2 discovered by van Velzen et al. (2010), for which 
VLA observations at v = 8.5 GHz place upper limits of Fy < 0.2 
mJy at ; ~ 7 and ~ 90 days after discovery. Unfortunately, these 
non-detections are not constraining because the flux is predicted to 
rise rapidly oc p prior to the peak at f ~ 1 year ( 34.2t . 

Radio emission might still be detectable now from candidates 
discovered in the past few years (e.g. Gezari et al. 2007, 2008, 
2009). Moreover, for a typical source redshift z ~ 0.2, the GHz 
emission from a tidal disruption could be sufficiently bright (~ mJy; 
eq. (9)) and extended that the diametric radio lobes (Fig.[T) physical 
(angular) separation ~ 2R„ ~ 0.5 pc (a 0. Imas), could be resolved 
using Very Long Baseline Interferometry, again analogous to the 
afterglow of GRB 030329 (Taylor et al. 2004). 

Another exciting prospect is the discovery of tidal disruptions 
with present or future radio surveys, such as Pi GHz Sky Survey 
(PiGS) of the Allen Telescope Array (Bower et al. 2010). Given 
PiGS predicted field of view (~ 10,000 deg^), cadence (~ days- 
years) and sensitivity limit ~ 5 mJy at v ~ 3.1 GHz, we estimate 
that ~ 100 tidal disruptions could be detected per year, assuming 
fiducial values for the jet properties from equation (|9]; a tidal dis- 
ruption rate ~ 10"^ yr"' per galaxy; and a local black hole density 
~ 3x 10"^ Mpc"' (e.g. Hopkins et al. 2007). Another possible strat- 
egy is a targeted survey towards nominally radio quiet spheroids or 
bulges, which are preselected to harbor black holes with masses 
Mbh ~ 10' - 10* Mo (D. Frail, private communication). 

6 DISCUSSION 

In analogy to other accreting black hole systems such as AGN and 
X-ray binaries, the tidal disruption of a star by a SMBH may be ac- 
companied by a powerful, relativistic jet. Such jet may be respon- 
sible to the acceleration of ultra-high energy cosmic rays (Farrar & 
Gruzinov 2009). In this Letter we have shown that the interaction 
of such a jet with the ambient ISM produces a bright radio-infrared 
transient that peaks on a timescale ~ 1 year after disruption. Ra- 
dio transients from tidal dismption may be detected blindly with 
upcoming radio transient surveys, or via follow-up observations of 
disruption events detected at other wavelengths on timescales of 
several months to years. Upcoming radio transient surveys similar 
to PiGS could discover tens or hundreds of tidal disruptions per 
year. The good radio localization would confirm the location of the 
transient at the center of the galaxy. 

The detection and characterization of tidal disruptions is a 
promising means to study the demographics of SMBHs (Gezari et 
al. 2009). Tidal dismption also provides a unique venue to study the 
physics of accretion, and its connection to jet formation, under con- 
ditions in which the 'boundary conditions' of the flow are relatively 
well-determined. The thermal disk emission provides a direct esti- 
mate of the accretion rate. Furthermore, because the magnetic field 
strengths of solar type stars are known, the magnetic flux through 
the accreting material is well constrained. 

An important unsolved question associated with jet formation 
is whether the magnetic field necessary to power the jet is advected 



with the flow (e.g. Spruit & Uzdensky 2005), or whether it is gen- 
erated locally in the disk by instabilities or dynamo action. Because 
the magnetic flux of a solar-type star is insufficient to drive a pow- 
erful jet, the detection of bright radio emission associated with tidal 
dismption would favor the hypothesis that locally-generated fields 
are responsible for jet (and vice versa). Furthermore, if only a frac- 
tion of (otherwise similar) tidal disruptions shows evidence for a 
jet, it would indicate that a second parameter in addition to the ac- 
cretion rate (such as the black hole spin) controls the jet strength. 
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